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Abstract. We discuss the limit of small width for the Laplacian defined on a waveguide with 
Robin boundary conditions. Under suitable hypothesis on the scaling of the curvature, we 
prove the convergence of the Robin Laplacian to the Laplacian on the corresponding graph. 
We show that the projections on each transverse mode generically give rise to decoupling 
conditions between the edges of the graph while exceptionally a coupling can occur. The 
non decoupling conditions are related to the existence of resonances at the thresholds of the 
continuum spectrum. 



1. Introduction 

The interest in the analysis of differential operators on metric graphs has been driven by the idea 
that some physical systems can be well modeled by using lower dimensional approximations. 
Intuitively one expects that graph-like approximations can be used to describe the dynamics 
in constrained systems characterized by two scales of length: a "large" one along the direction 
of the edges and "small" ones in the transverse directions. 

In mathematics a metric graph is a one dimensional singular manifold, and is defined by as- 
signing points, the vertices connected by a set of oriented segments, the edges. The dynamics 
on the graphs is fixed by defining a differential (or pseudo-differential) operator on the graph 
and boundary conditions in the vertices. In the following we shall consider only self-adjoint 
operators. 

When this kind of structures arise as approximations of quantum systems, it is customary to 
call them quantum graphs. In this case the self-adjointness assumption is natural. One of the 
first and most famous applications of quantum graphs dates back to 1953, when they were used 
to model the dynamics of 7r electrons in organic molecules (37] . More recently a renewed interest 
in quantum graphs has resulted from the development of nanotechnologies. At present devices 
based on carbon nanotubes and metallic nanowires are commonly produced and studied. In 
such structures the mean free path can reach hundreds of micrometers , while the transverse 
confinement can be of the order of ten nanometers. Because of their small dimensions and 
purity such devices represent an ideal framework for the analysis of many peculiar phenomena 
of quantum mechanics, see, e.g., [IS] and [29], and in many applications they can be treated 
as one dimensional systems. For a comprehensive review on properties and applications of 
quantum graphs we refer to [21] , [25] , [26] and [1] . Even in the field of classical mechanics there 
is a huge number of problems in which metric graphs define simplified but non trivial models. 
Typical examples arise in the analysis of electromagnetic or acoustic waves in thin waveguides. 
Many efforts have been done in the last fifteen years to understand to which extent a one 
dimensional dynamics on a metric graph approximates the dynamics in a constrained system. 
Aim of this work is to investigate some relevant features of graph-like approximations for the 
quantum dynamics in networks of thin tubes. 

To understand the core of the problem it is sufficient to discuss the case of a graph with one 
vertex and iV infinite edges. The natural Hilbert space for such a system is the direct sum of 
iV copies of L 2 ((0, oo)) and its generic element is (fi, . . . , /at) with e L 2 ((0, oo)). We restrict 
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ourselves to a setting in which the dynamics on the graph is generated by an Hamiltonian that 
on each edge coincides with the (positive) Laplacian. From the mathematical point of view 
one can define several self-adjoint operators on the graph that coincide with the Laplacian on 
the edges. Each element of this family of operators is identified by the boundary conditions in 
the vertex imposed on the functions in its domain. In the following we shall call the boundary 
conditions in the vertex gluing conditions. A linear relation between (/i(0), . . . , /tv(0)) and 
(f[(0), . . . , f' N (0)) must be used to fix the gluing conditions in a way such that the corresponding 
operator on the graph is self-adjoint. Making use of Krein's theory one can characterize all the 
possible self-adjoint gluing conditions in the vertex. This was done by Kostrykin and Schrader 
in [21] (see also [16]); for each vertex of degree N there are N 2 real free parameters to fix the 
gluing conditions. 

Two well known examples of self-adjoint gluing conditions are 
(1.1) /;(0) = Vj = l,...,iV 

and 

N 

(1-2) A(0) = / a (0) = • • • = MO) , # (°) = • 

3=1 



Condition (1.1) is usually called decoupling condition or Dirichlet condition. We shall use the 
expression decoupling condition to keep in mind that this type of gluing in the vertex implies 
decoupling among the edges, i.e., no transmission trough the vertex is possible. Condition (1.2) 
is usually referred to as condition of free type or Kirchhoff type. In our opinion the expression 
condition of Kirchhoff type is a bit misleading. It recalls current conservation in electric circuits 
but as a matter of fact every self-adjoint gluing condition conserves the quantum probability 
current across the vertex. The expression free condition seems more appropriate because (1.2) 
generalizes the free one dimensional Laplacian to a non trivial topology. 

In view of applications one is interested in understanding which gluing conditions are more 
appropriate to model strongly constrained quantum systems. A standard strategy to approach 
this problem consists in the analysis of the limit, in some suitable sense, of the Laplacian in a 
network of thin tubes as the network shrinks to the underlying graph. 

Let us assume that the network we consider is made up of tubes that far from the vertex are 
straight and of constant width. To fix ideas we also suppose that the tubes have all the same 
width. In this setting, far from the vertex, the dynamics is factorized in the direction along the 
axes of the tube and in the transverse direction. Moreover as the manifold shrinks the energy 
gap between the transverse modes increases as the inverse of the squared width of the tube. 
Even if this simple picture fails as one approaches to the vertex, it suggests that the natural 
way to reduce the dynamics to one on the underlying graph is to project onto the transverse 
modes and that, in the limit of zero width, each projection can lead to a unitary dynamics, 
that is an effective dynamics which leaves invariant the subspace associated to the transverse 
mode. 

The feasibility of this procedure and the corresponding limit operator on the graph depend 
on the boundary conditions that one imposes on the boundary of the tubes and on which 
transverse mode the projection is taken. 

In this paper we discuss this problem in the most simple geometrical setting. We take as 
initial domain Q, a strip of constant width d around a base curve T and we assume that T has 
no self-intersections. Under this assumptions the model is greatly simplified by the presence 
of a global system of coordinates (s, u) adapted to the geometry of the system, where s is 
the arc length coordinate on V and u is the orthogonal coordinate. The relevant geometric 
quantity are d and the scalar curvature of the curve, 7(s), that in our hypothesis is a smooth 
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and compactly supported function. We rescale the initial domain according to d — > e a d and 
7(s) — > e~ lr y(s/e), where e > and a is a big enough positive constant, and we obtain a 
net of domains Q £ which as e — > collapses onto a broken line, i.e., onto a graph made up 
of one vertex and two infinite edges. We consider the Laplacian on Q £ with Robin boundary 
conditions and discuss its convergence to the Laplacian on the limit domain which can be seen 
as the most simple example of graph. 

The most relevant technical part of the present paper consists in the proof that the projection 
of the dynamics on the n-th transverse mode is unitary in the limit e — > 0, and that the effective 
dynamics is described by the Hamiltonian 



d 2 6 
h £ = -— + —>y 2 (-le) 
ds 2 + e 2lUh 



where (3 n are some coefficients related to the energy of the transverse mode. 
The analysis of the limit for e — > of Hamiltonians of the form of h e n was performed in the 
former work by the same authors [3 J . For each n the limit operator depends on the low energy 
properties of the Hamiltonian h n 

d 2 

= 'ds 2 + ' 

More precisely, under our assumptions on 7(5) two cases can occur 

(1) There does not exist a zero energy resonance^] for h n , then h e n converges to the Laplacian 
on the graph with decoupling gluing conditions in the vertex. 

(2) There exists a zero energy resonance f r>n , for h n . In such a case one can define two real 
constants c± = lim^-too f r , such that c 2 + + c 2 _ = 1, and the limit operator on the graph 
is the Laplacian with gluing conditions given by 

(1.3) c_A(0) = c+/ 2 (0) , c + /((0) + c_£(0) = . 

The convergence has to be intended in the norm resolvent sense. 

We remark that the existence of a zero energy resonance is an exceptional event, that in our 
case is related to some very special choices of the initial curve T, see [3] for few examples. This 
implies that in most of the cases the limit operator is defined by decoupling conditions in the 
vertex. For this reason we call the case 1 generic and the case 2 non-generic. 



The gluing conditions (1.3) are known in literature as scale invariant, see [T7j, and are param- 
eterized by one independent real parameter, e.g., the ratio c + /c_. By using the result proved 
in [7] we shall see that, in the non-generic case, a deformation of order e of the angle 9 between 
the edges of the graph leads to a more general class of gluing conditions, i.e., the ones defined 
by 

c_A(0) = c + / 2 (0) , c + f[(0) + c_£(0) = 6(c+A(0) + c_/ 2 (0)) 
where b is a real constant related to the deformation parameter. 

It is worth noticing that in the non-generic case the gluing conditions in the vertex implies 
a coupling between the edges; our analysis includes a wide class of boundary conditions and 
holds for all transverse modes. 

The first results in the same direction presented here, in a setting with several edges, concerned 
the case of compact networks of tube^Jwith Neumann boundary conditions. In [13J was proved 
the convergence of the solutions of the heat equation on the network to the corresponding 
solution of the equation on the graph. The proof made use of some results on the convergence 



1 See section |3] for the definition of zero energy resonance. 

2 For compact network we mean a network that is contained in a compact region of the space. For such a 
network and for the corresponding underlying graph the spectrum of the Laplacian is discrete. 
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of Markov processes proved in the same paper. In a similar setting the convergence of the 
spectrum was proved in [27], [36], [38], [2H] and [TT]. In all these works the gluing conditions 
arising in the limit are of free type. The most recent result on the Neumann problem (and 
in a setting in which the approximating manifold has no boundary), was given in [34]. In the 
latter work in the case of compact and non compact networks, was proved the strong resolvent 
convergence to the operator on the graph with gluing conditions in the vertex of free type. In all 
the aforementioned papers, only the projection onto the lowest transverse mode was considered. 
The case we just described is the most simple one. In this setting the energy of the lowest 
transverse mode is equal to zero and, in a neighborhood of the vertex, one is allowed to ap- 
proximate the wave function on the network with a constant function. The case with Dirichlet 
boundary conditions has revealed much more tricky. This is due to the fact that in the latter 
case the energy of the transverse modes always increases as the inverse squared width of the 
tubes and this forces to rescale the Hamiltonian by subtracting the divergent energy term. The 
scaled operator has a finite number of eigenvalues which all diverge as the network squeezes. 
A first result on the problem with Dirichlet boundary was given in [33] . In a setting of a compact 
network it was proved that the limit operator on the graph is characterized by decoupling 
conditions in the vertex. It is worth noticing that in this work the decoupling was obtained as 
a consequence of an ad hoc hypothesis on the volume of the manifold in a neighborhood of the 
vertex. 

Recently D. Grieser [15] has proved that, for a large class of boundary conditions, generically 
the limit gluing conditions are of decoupling type. This was already argued by S. Molchanov 
and B. Vainberg, see [30]. In these works, for compact networks, the spectral convergence of 
the Laplacian on the network to the operator on the graph is studied; the approach is based 
on the analysis of the scattering problem associated to the network of tubes and makes use of 
the analytic properties of the resolvent of the Laplacian on the manifold. The special cases 
in which the coupling occurs are related to the existence of singularity of the resolvent at the 
thresholds of the energy of the transverse modes. 

Earlier the existence of a non-decoupling limit in the Dirichlet case was proved in [3] for the 
same model discussed in this paper. Such a model was proposed for the first time as a prototype 
of a Dirichlet network collapsing onto a graph in [8] where the generic case leading to decoupling 
conditions was discussed. 

In view of applications the most relevant type of boundary conditions in the modeling of con- 
strained quantum mechanical systems is the Dirichlet one while Neumann boundary conditions 
arise mostly in the case of electromagnetic or acoustic waveguides. Robin conditions are used 
for example in numerical simulations to model the interface between semiconductors, see, e.g., 
[39]. 

The interplay between geometry and boundary conditions in waveguides has been studied in 
many works (see, e.g., [12J, [10], [9], [6], [14] and references therein). Most of them focus on the 
differences between Neumann and Dirichlet boundary. By changing the Robin constant one 
can continuously switch from Neumann to Dirichlet boundary conditions. For this reason we 
guess that the analysis carried on in this paper can help to gain a deeper understanding in this 
problem. 

Let us stress that we consider only "symmetric" waveguides, i.e., we take the same boundary 
conditions on the upper and lower boundary of the waveguide. Asymmetric boundary conditions 
would lead to a divergent term in the transverse Hamiltonian of the order of the inverse of the 
width of the waveguide. Some results on the spectral properties of the Laplacian in asymmetric 
waveguides can be found in [22] and [23] . 

The present paper is structured as follows. In section [2] we define the Laplacian with symmetric 
Robin boundary conditions on the waveguide and introduce the correct scaling to get a family 
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of waveguides that collapses onto a graph. In section [3] as preliminary results, we discuss 
the spectral structure of the one dimensional Laplacian on a compact interval with Robin 
boundary conditions and we recall some results on the limit of Hamiltonians with short range 
scaled potentials taken from [3]. After that we state the main theorem. Section [4] is devoted to 
the proof of the main theorem. In section [5] we discuss the effect of small deformations of the 
relevant parameters of the problem. In this section we make use of the result proved in [7j where 
the same problem was studied in the case of a Dirichlet boundary. A section of conclusions 
and remarks closes the paper. The proofs of few technical estimates and of a resolvent formula 
that will be used to prove the main theorem are postponed in appendix |A| 

2. The model 

Let T be a curve in IR 2 given in parametric form by T := {(71(5), 72(5)), s G IR} and let us 
assume that it is parameterized by the arc length s, i.e. 7^(s) 2 + 7 2 (s) 2 = 1. The curve T is 
completely defined up to isometries once the signed curvature 7 is known 

the curvature radius of T in s is equal to the inverse of the modulus of the signed curvature. 
We shall assume that 7(5) G Cq°(R), therefore T is a straight line outside a compact region. 
We shall also assume that T has no self-intersections. Thus T consists of two straight lines, li 
and I2, with the origins, 0\ and O2, connected by an infinitely smooth, non self-intersecting, 
curve C, running in a compact region. The integral of 7 gives the angle 9 between l\ and I2 

(2.1) 9 = [ «f(s)ds. 

Jr 

Let us denote the open strip of width 2d around T by Q: 

Q := {(x, y) s.t. x = 7i(s) - wy' 2 (s), y = 7 2 (s) + wy[(s), s G IR, u G (-d, d)} . 

We assume sup s \j(s) \d < 1, in this way (s,u) provide a global system of coordinates in Q. 
Let us define the sesquilinear form on L 2 (Q) x L 2 (Q) with domain 

given by 

Q n [if, ip] := dx dyVcp Vtp . 
Jn 

It is well known that Qn is closed and positive and that the associated self-adjoint operator is 
the Laplacian in the domain Q with Neumann boundary conditions. 
Now we consider the following perturbation of Qq depending on a G IR 

:= / dxdyW^ViP + a I (^^)\ m dS 
Jn Jan 

where dS is the Lebesgue induced measure on dfl. Then by Sobolev embedding theorems (see, 
e.g., [T]), is a small perturbation of Qq in the sense of quadratic forms and is closed 
and bounded from below on 

@(Qq) := H\Q) x H\Q) . 
Moreover C7 CO (M 2 ) x C^°(1R 2 ) is a core for Q n and (see, e.g., [35j Th. X.17). 
We denote by — the self-adjoint operator associated to Q^. One can verify that the operator 
— Aq coincides with the Laplacian with Robin boundary conditions on dQ, i.e., functions in 
Aq) belong to H 2 (Q) and their trace on dfl satisfies the boundary condition ^\an + 
oti>\an = 0. 
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ip e H 2 (tt) s.t. 



dip 
dn 



+ aip 



mi 



mi 



and let us recall that by the first representation theorem of quadratic forms, see [20], we have 

#(-A£) := {iP G ®(Q%) s.t. 3 X G L 2 (fi),Vy? G #(fig) , Qgfo ^ = (^x)^(0)} • 

Integrating by parts we have immediately that & C @(— A^). Now we prove the reverse 
inclusion and then equality follows. Let us assume that ip G £^(— Aq) then there exists x G 
L 2 (f2) such that 

Qnb.V'] = 0,x)i2(fi) 
for all (y? G Co°(n). For such a ip we simply have 

yRr 



q^pM = / dxd y v<pVip 



and then 



dx dy (px 



dx dyV(f Vip 



dx dy ip (—Aip) 



which implies that — Aip G L 2 (VL) that is ip G H 2 (£l). Now we take ip G ^(Qq) and using 
Gauss-Green theorem, we find 



GSM 



(ix (iy ( — A^) + / cLS 1 </j| 

9f2 



r«2 



90 



9n 



+ aw 



mi 



mi 



This is a bounded functional with respect the L 2 (Q) topology in (p if and only if the boundary 
conditions ^\ 9n + aip\ dn = hold and then 2>(-A^) = 

In order to study the properties of — A^ it is convenient to use the coordinates (s,u) which 
belong to Q' = R x (—d, d). The following proposition holds true. 

Proposition 1. For 7 G Cq°(R), — Aq is unitarily equivalent to the operator H in L 2 (Q', ds du) 
defined by 
(2.3) 



dip 



dip 



ip G H 2 (Vt') s.t. — -(s, d) + a x (s)ip(s, d) = , —£-(s, -d) + a 2 (s)ip(s, -d) = 



where 
(2.4) 

and 

with 
(2.5) 



9m 



a — 



9m 



7(s) 



if :-- 



2(1 + dj(s)) 




a 2 (s 



a + 



7(s) 



2(1 -d 7 (*)) 



1 



9 9 2 



9s (1 + m 7 (s)) 2 9s 9m 2 
7 (s) 2 ui'(s) 5 M 2 7 '(s) i 



+ 



4(l+w 7 (s)) 2 2(1 + m 7 (s)) 3 4(l + w 7 (s)) 4 ' 
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Proof. We denote by Qq the sesquilinear form obtained by restricting Qq to C7^°(IR 2 ) x l7^°(IR 2 ); 
the closure of Qq is Qq. The Hilbert space L 2 (Q,dx dy) is mapped by the change of variables 
into L 2 (Q', (1 + wy)dsdu). With these coordinates the form Qq reads 

+ a [(l + dj(s))lp(s,d)il)(s,d) + (l — d / y(s))'(p(s,—d)ip(s,—d))jds. 



Notice that even in the new coordinates the domain of Qq is C^(IR 2 ) x C^(]R 2 ). We consider 
the unitary map U : L 2 (Q', (1 + wj)ds du) —>■ L 2 (Q', ds du) given by 

(2.6) (Uij)(s,u) := (l + u-f(s)) 1/2 ij(s,u). 

A straightforward calculation shows that the form Qq is unitarily equivalent to the form Qq, 
in L 2 (Q', ds du) defined by 

9{Q&) := C™(R 2 ) x C °°(R 2 ) , 
+ / (a 1 (s)'<f(s,d)ip(s,d) + a 2 (s)'<f(s,—d)ip(s,—d))ds 



where «i(s) and a 2 (s) are given by equation (2.4) and K by equation (2.5). Now we just need 
to compute Qq,, the closure of Qq,. Let us consider 



Qcv[<P,if>] 



1 



d(p dip d(pdip\ , 
■ "t; — -7; — I- — t; — ds du 



[1 + wy) 2 ds ds du du J 

with the same domain as Qq,. Since < c < (1 + dj) < c" 1 for some positive constant c, we 
notice that QQi[ip,ip] is equivalent to the H 1 norm in Since V and «j are bounded, then 
Qq, is a small perturbation of Qqi and its closure is given by 



+ / (ai(s)ip(s,d)ip(s,d)+a2(s)ip(s,—d)ip(s,—d))ds. 
Jr 

It is easy to repeat the argument used before to prove that the domain of — is equal to 3> 
defined in (2.2), to see that the self-adjoint operator associated with Qq, coincides with the 
operator H given by (2.3) - ( |2.5 ). □ 

From now on we denote L 2 (Q' ,ds du) simply by L 2 . 

In order to get a family of waveguides that collapses onto a prototypical graph, we rescale the 
geometric parameters of the system 7 and d in the following way: 

(2.7) j(s) — - £l (£) 

(2.8) d — ► 5 £ d. 

Where e > and 5 £ > are two dimensionless scaling parameters, such that 5 s — > as e — > 0. 
With this scaling we have two characteristic scales: e is the scale of variation of the curvature 
7 while 5 e is the intrinsic scale in the transverse direction. We shall consider the regime where 
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5 £ <C e, that is, the curvature is slowly varying with respect to the width of the waveguide. In 
particular we assume 5 £ := e a with a ^ 1 big enough, a more precise statement on the ratio 
5 £ /e will be given in the following. Notice that the angle 9 defined in (2.1 ) is unchanged by the 
scaling (2.7). 

We also rescale the Robin constant a in the following way 

e a 
a — > a := — , 

5 £ 



in this way the Robin boundary conditions in (2.2) are invariant under the scaling of d. 
We obtain a family of domains fl £ and of operators — such that Q £ approximates, for 
e —>■ 0, the broken line of angle 9 made up of two half lines, l\ and I2, with the same origin, 
0\ = O2 = O. Since we assume a ^ 1 then (s, u) are a system of global coordinates also for 
Q £ . 

From proposition [l] it follows that, for all e > 0, the operator — is unitarily equivalent to 
the self-adjoint operator H £ : S#(H £ ) C L 2 — > L 2 given by 

9 1 9 1 d 2 1 



with 



;= 7 (^) 2 + 8 £ /eui'{s/e) 5 (F/e) V 7 '(s/e) 2 



4(l + urf(s)) 2 2(l+mf(s)) 3 4 (l + wrf(s)) 4 

and with domain 
(2.10) 



^(# £ ) := |^ e # ) s.t. ^-(s, d) + af (s)V»(s, d) = , — ^-(s, -d) + al(s)V»(s, -d) = 
here af(s) and a|( s ) are gi ven by 

(2.11) of (a) :=a- -r, ^ . .. , a|(s) := a + 



2(1 + drf(s)) ' ^ w ' 2(1 -drj £ (s)) 

and we have introduced r] £ (s) := 5 £ /e , ~f(s/e). In the following 77 s will play the role of a small 
quantity in a suitable topology and will allow a perturbative analysis. 



3. Main result 



Equation (2.9) shows that the transverse kinetic energy is divergent in the limit e — > 0. This 
is a common problem for these kind of singular limits. In order to overcome this problem, it 
is convenient to introduce an s dependent orthonormal complete set of states of L 2 ((—d,d)) 
which diagonalizes the transverse part of the Hamiltonian and provides a useful framework to 
discuss the limit of H £ in the sense roughly described in the introduction. To this aim we start 
this section with a short discussion on the one dimensional Robin Laplacian in L 2 ((—d, d)). 
Given two real constants ot\ and a>2, we denote by h ai<a2 the Robin Laplacian on L 2 ((— d, d)); 
h Qll a2 is the self-adjoint operator defined as 

@(h ai>a2 ) := G H 2 ((-d, d)) s.t. V'(d) + a^{d) = , -ip'(-d) + a 2 ip(-d) = 0} 

d 2 if) 



GRAPH-LIKE MODELS FOR THIN WAVEGUIDES WITH ROBIN BOUNDARY CONDITIONS 9 

Let us denote by g ai ,a 2 (k 2 ) '■= (hai,a 2 — k 2 ) -1 the resolvent of h QljQ , 2 . The integral kernel of 
gai,a 2 (k 2 ) is explicitly known: 

1 sin[k(2d - \u - u'\)\ 



9a±,a 2 {k 



(3-1) 



2kd? cos{2kd) 

k(ai — a 2 ) sin[k(u + u')} — (aia 2 + k 2 ) cos[k{u + u')} 
2kd 2 [(a\a 2 — k 2 ) sin(2kd) + k(a± + a 2 ) cos(2kd)} 
(a\a 2 — k 2 ) cos[A;(u — u')] 



2kd 2 cos(2A;(i)[(a;ia2 — k 2 ) sin(2kd) + k(ai + a 2 ) cos(2kd)] 
for k 2 G p(h ai ^ a2 ) and Im k ^ where p(h Ql]Q2 ) denotes the resolvent set of h QliCf2 . 

We denote by A n , n — 0, 1, 2, . . ., the eigenvalues of h aiiQ2 arranged in increasing order. Using 

(3.1) it is straightforward to prove that X n = k 2 with k n given by the solutions of 

(3.2) («i«2 — k 2 ) sm{2k n d) + k n (a\ + a 2 ) cos(2k n d) = 0, 

positive eigenvalues correspond to k n e IR + , while negative eigenvalues are given by k n e iM. + . 
The corresponding eigenfunctions have the form 

(3.3) 4>n(u) = A n sin(fc n w) + B n cos(k n u) n = 0, 1, 2, . . . , 

where A n and B n are suitable coefficients. The eigenfunctions <p n can be chosen real. When we 
want to stress the dependence on the boundary conditions, we shall denote the eigenvectors of 
hai,a 2 by 0n(a, u) , where a = (ai,a 2 ). 

We simply denote h a ^ a by h a and by /i n , n — 0, 1, 2, . . ., its eigenvalues arranged in increasing 
order. The eigenvalues of h a can be written with p n satisfying: 

(3.4) p n sin(p n d) — a cos(p n d) = n — 0,2,4,... 

(3.5) p n cos(p n d) + a sm(p n d) = n = 1,3,5, ... . 

Positive eigenvalues correspond to p n G M + , while negative eigenvalues are given by p n G 
The corresponding eigenvectors have now a definite parity and can be written as 

(3.6) £ n (u) = N n cos(p n ii) 7i = 0,2,4,... 

(3.7) £ n («) = N n sin(p n u) n = 1,3,5,..., 

where N n is the normalization constant. For a ^ all the eigenvalues are non negative, 
for — 1 ^ ad < there is one negative eigenvalue and for ad < — 1 there are two negative 
eigenvalues. In figure [l] the first four eigenvalues of h a are plotted as functions of a (for d = 1). 

Let us now take «i and a 2 in h Qlia2 of the following form: 

i] i] 
a± := a — — — , a 2 := a + 



2(1 + drj) ' 2(1 - dr?) ' 

For rj <C 1 the Hamiltonian h Ql Q , 2 can be considered as a perturbation of h a . The eigenvalues 
of h aia2 , X n = k 2 , are defined by the equation A(k n ,r)) = with 



a + 2Ti^)( a -2Tr^) 1 - fc 



+ k n 



2(1 -drj) ) V 2( 1 + ^)/. 



sin(2A; n (i) 
cos(2A; n <f) . 
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Figure 1. Plot of the first four eigenvalues of h a as functions of a, in the plot 
it is assumed d = 1. 



We want to give a perturbative expansion of the eigenvalues A n in the small parameter 77 up to 
the second order: 

2 



3=0 

The length d appears in the expansion for dimensional reasons only, the small parameter here 
is rj. The coefficients kn , for j = 0, 1, 2, in the expansion k n = k { ° ] + k^d V + k ( n\dri) 2 + Q{rf) 
can be obtained from the expansion A(k n ,r)) = A^ + A^dr] + A^ 2 \dr]) 2 + 0(r] 3 ) and imposing 
A^' = 0, j = 0, 1, 2. A straightforward calculations gives 



1.(0) 
1.(1) 

1.(2) 



Pr, 




p n [a-2d(a 2 +p 2 )] 



4d 2 (a 2 +p 2 n ) [a + d(a 2 + p 2 n )] 
where p n were defined in equations (3.4) and ( |3.5[ ), and since A n = k 2 we immediately obtain 

A« = 

li n [a - 2d(a 2 + /x n )] 



(3.8) 



A[ 2 > 



The coefficients An 



(2) 



n 



0,1, 



2rf 2 (a 2 + y u n )[a + rf(a 2 + y u n )] ' 
., are implicitly defined as functions of a. 



Equations (2.9), (2.10) and (2.11 ) indicate that for each s the transverse kinetic term of H £ has 



the same form of h aijQ , 2 and that our hypothesis allows a perturbative analysis of its spectrum. 
Therefore for all s e R we consider the Hamiltonian h a e a e ; L 2 ((—d,d)) — > L 2 ((—d,d)), 



where af(s) and a|(s) were defined in (2.11). Such Hamiltonian depends on £ and s only via 



the constants af(s) and in the boundary conditions. In the following we shall use the 



notation a £ (s) = 



»)■ 
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The eigenvalues and the eigenvectors of h a ^ a e are defined according to ( 3.2[ ) and (3.3) and will 



be denoted by A^(s) and <^(s). We shall use also the notation <p £ n (s) = </> n (a e (s)) to remind 
that <f>n(s) depend on s and e only via a\(s) and a|(s). For fixed sel, {0n(« e (s))} n _ o 1 is 
an orthonormal basis of L 2 ((—d,d)). 

For all z G C\R and for all m, n = 0, 1, 2, . . . , we denote by r £ mn {z) the reduced resolvent of 
H £ , i.e., the operator in £g(L 2 (R), L 2 (R)) defined via its integral kernel b}| 

r £ mn (z]s,s) := I I (f) m (a £ (s),u) [h £ - - z^j (s,u; s' ,u) <p n (a £ (s'),u)dudu . 



-dJ-d 

e2 



Notice that we have subtracted the divergent quantity fi n /5 £ from H £ in order to compensate 
the divergence of the transverse kinetic energy and get a non trivial limit; this procedure was 
already used in [33], [8], [3] and [?]. 



We also need to recall some facts on one dimensional Schrodinger operators with short range 
potentials in L 2 (IR). Let us consider the Hamiltonian h given by: 

(3.9) h: =-at + v{s) > 

and let us assume that for some c > 



(3.10) / v(s)ds ^ e c{ - l v G L x { 



We say that h has a zero energy resonance if there exists f r G L°°(R), f r ^ L 2 (IR) such 
that hf r = in distributional sense. Furthermore, if f r exists, it is unique, up to a trivial 
multiplicative constant and one can define two constants 

(3.11) c_ := lim f r (s) and c + := lim f r (s) . 

The constants c_ and c + can not be both zero, in such a case f r would be in L 2 (IR), then 
zero would be an eigenvalue for h (see Lemma 2.2. in [5]), but this is impossible under our 
assumptions on v, see Theorem 5.2. in [19]. We can choose c_ and c + real and such that 
c 2 +c 2 =l. 

Let h r be the following family of self-adjoint operators depending on c_ and c + 

(3.12) ®(hr) := {/ EH 2 (R\ 0) s.t. c„/(0 + ) = c + f(0~) , c + f'(0 + ) - c^f(0~) = 0} 

(3.13) hr f,= - d l^ s ^o. 

The Hamiltonian h r is a self-adjoint extension of the symmetric operator —A in dimension one 
defined on Cq°(R \ {0}). For c_ = c + the operator h r coincides with the free Laplacian on 
the line; we refer to [2] for a comprehensive characterization of the point perturbations of the 
Laplacian in dimension one. Let us notice that the operator h r can be rewritten as an operator 
on if 2 ((0,oo))©# 2 ((0,oo)) C L 2 ((0,oo))©L 2 ((0,oo)) by defining f x and f 2 in F 2 ((0,oo)) such 
that fi(x) = f(x) for x > and f 2 (—x) = f(x) for x < 0. Within this notation the condition 



in x = in the domain of h r can be written as we did in the introduction, see equation (1.3). 
We denote the one dimensional Laplacian with decoupling (or Dirichlet) gluing conditions in 
the origin by ho 

@(h ) := {/ G H 2 (R \ 0) D H^R) s.t. f(0) = 0} 



3 We denote by S§(L 2 {J-L),L 2 (TL')) the Banach space of bounded operators from Ti — > H! and by 
\\3g{L2(H),L?CH')) the corresponding norm. 
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d 2 f 



hof 



ds 2 



Even the operator ho can be written in the standard notation of the Laplacian on graphs, the 
condition in x = would correspond to a decoupling condition in the vertex, see equation (1.1 ). 
Now we rescale h in the following way 



(3.14) 



d 2 1 / i \ 



and we discuss the convergence of h £ in resolvent sense. The following proposition is taken from 
lemma 1 in [3]. 



Proposition 2. Take h and h £ defined as above and assume (3.10). Then two cases can occur: 
(1) There does not exist a zero energy resonance for the Hamiltonian h, then 



u — \\m{h £ — z) 1 = (ho — 



£^0 



z e C\ 



(2) There exists a zero energy resonance f r for the Hamiltonian h, then 



u — lim(/i £ 

£->0 



(h r 



zeC\ 



where h r was defined in (3.12) and (3.13). 
We denote by (3 n the coefficients 

(3.15) (3 n :=-l/A + \^ 

and by h n the Hamiltonian 



h n :-- 



'ds 2 



+ M 2 {s). 



Our main result is stated in the following theorem. 



Theorem 1. Assume that V has no self-intersections and that 7 G C ° 
a > 3, then for all n, m — 0, 1, 2, . . . two cases can occur: 

(1) For all n such that there does not exist a zero energy resonance for h n we have 

-1 



moreover take 



u — lim r\ 

e— >0 



O~m,n{ho — Z) 



zeC\ 



(2) For all n such that there exists a zero energy resonance, f Tyn , for h n we have 

z eC\R. 



u 



\imr £ mn (z) = S mn (h rtn - z) 1 



where h rn is defined according to equations (3.12) and (3.13). 

Notice that the constants c_ )n and c +jn in the definition of £^(/i rjra ) may depend on n and are 
related to f rn via equation (13.11 ). 



We remark that the parameter a in the boundary conditions affects the limit only through the 
coefficients (3 n (a). In figure [2] the functions (3 n (a) are plotted for n = 0,1,2, 3 and d = 1. The 
exceptional case f3 n (a) = is included in the statement of the theorem because in such a case 
the Hamiltonian h n has a zero energy resonance given by the constant function, therefore c_ )n 
and c +j „ coincide and the limit operator h r ^ n is the free, one dimensional, Laplacian. 
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Figure 2. Plot of fiJa) for n — 0, 1, 2, 3 and d= I. 



4. Proof of theorem \T\ 

Before getting into the technical core of this section devoted to the the proof of the main 
theorem, let us spend few words on the strategy we shall follow. In our regime the curvature 
is slowly varying with respect to the transverse dynamics. In particular this means that cnf(s) 
and a|(s) are slowly varying functions and that each subspace corresponding to an eigenstate 
(f) e n is "adiabatically protected" in the limit e — > 0. This allows us to split the proof into two 
steps. 

(1) First we prove that each subspace corresponding to an eigenstate 4> £ n is adiabatically pro- 
tected in the limit e —>■ and we prove that the leading term in the reduced Hamiltonian 



(up to the renormalization factor ji n /8 e ) is 



(4.1) 



hi 



ds 2 



i @ n 2/ / \ 

+ -,7 ( s / £ ) 



where (3 n were defined in (3.15). This is done in lemma [T] and lemma [2] 
(2) As a second step we study the limit of h e n . Here we shall make use of the proposition 
[2] to prove that, for each n, two cases can occur: if the potential /3„7 2 generates a zero 
energy resonance for h n , then h £ n converges to an operator of the family defined in (3.12) 
- ( |3. 13 ) otherwise the limit operator is h , i.e., the Laplacian on the line with decoupling 
conditions in the origin. 

In the proofs c will denote a generic positive constant whose value can change from line to line. 
Let H e be the Hamiltonian 
(4.2) 



dip 



dip 



ip e H 2 (Q') s.t. -4s, d) + al(s)ip(s, d) = , —£-(s, -d) + a £ 2 {s)ip{s, -d) = 



du 



du 



(4.3) 



H £ 



2 



1 d 2 



1 7 2 ( S /e) 



ds 2 5 2e du 2 
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For all z e C\R and for all m, n = 0, 1, 2, . . . , we denote by r £ mn {z) the reduced resolvent of 
H £ , i.e., the operator in =^(L 2 (IR), L 2 (IR)) defined via its integral kernel by 

r-d 




dJ-d 



>„,{<r i >> ). '/ ) { H £ — — z ) (s, w; s', w') (p n (a £ (s'),u')du du . 



In the following lemma we prove that f^ n approximates r e mn . 
Lemma 1. Lei 7 G Cq°(R) and a > 3 then for all n, m — 0, 1, 2, . . . 

u - lim (r^ n (z) - f^ n (z)) =0 V.GC\i. 

Proof. It is sufficient to prove that there exists e such that for all < e < e and for all 
f,geC?{R) 

(4-4) IfeKnW-^nW]/), 



^ ce a 3 nc/i| L 2 



We shall make use of the fact that given a self-adjoint operator A in some Hilbert space 7i the 
inequality 

(4-5) ' 



\\(A-z)- 1 \\ a( H 7{) < 



Im 2 



holds. From the first resolvent identity and from the definition of r £ mn {z) and r £ mn {z) we have 
that 

-1 d 2 / . 



+ 
+ 



H £ 





-) 






-') 




6 e2 


-') 



<%( ff -^--~ 



-1 1 







rv 










/ L 2 



-1 



L 2 



-1 



L 2 



where 



W £ (s,u) 



b\(s,u) ~- 
6 £ 



_2u(5 £ /e)-f(s/e) + u 2 (5 £ /e) 2 -f(s/e) 
(l + ur] £ (s)) 2 

-l(s/e) 2 2 ul (s/e) + u 2 (5 £ /e)"t(s/e) 2 



b £ 2 (s,u) 



2u(5 £ /e 2 )i(s/e) 



+ 



(l + UT] £ (s)) 3 ' 

ui'{s/e) hu 2 (5 £ /e)i{s/e) 



4 (1 + U7] £ (s)) 2 2(1 +ur] £ (s)) 3 4 (1 + wrf(s)) 4 

6f , 6| and W e are bounded functions; more precisely, there exists eo such that for < e < eo 

5 £ 



(4.6) 

We use the notation 



< c- 



II^IU- <c^; H^IUoo ^c- 



i? e (^) := (H £ -z)~ l 



From the Cauchy-Schwarz inequality and estimates (|4.5|) and (4.6) we have 

d 2 



(4-7) |(<7,[r^)- f m>)]/)* 



Im z\ 



e 



ds 2 



R £ [z+ r " 



S s2 



L 2 
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The estimate of the third term on the right hand side of (|4.7|) comes directly from (|4.5|) 



(4.8) 



R £ z + 



$e2 



< \lm z\ 1 



L 2 



In order to estimate the other two terms we prove that there exists s$ such that for < e < Eq 

52 , ,, \ -1 



(4.9) 



ds 2 \ 5 e2 



z 



sC ce 



L 2 



To this aim we make use of an explicit formula for the resolvent R £ (z). The proof of this 
formula is postponed to appendix |A| 

We denote by H £ the self-adjoint operator in L 2 with the same formal expression as H £ but 
with domain characterized by boundary conditions not depending on the s variable 



(4.10) @(H £ ) := \ i) G H 2 (Q')s.t.-^(s,d) + a^(s,d) =0 



du 



du 



(s, —d) + aip(s, —d) = > , 



(4.11) 



H £ :- 



d 2 



l d 2 



1 ^{s/e) 



ds 2 5 2e du 2 



Moreover for all z G C\R we define 



R&z) := (H e - zY l . 

It will be crucial in the following that the derivatives of Rq{z) with respect to s and with respect 
to u commute. In particular we notice that — — ^ 7 and — prj^ commute with Rq(z). 
Let us introduce some notation and state few preliminary results. Let L 2 (dQ') = L 2 (IR) ©L 2 (1R) 
and q G L 2 (IR) © L 2 (M) denote a couple of functions qi G L 2 (R) for i = 1,2, q has to be 
understood column vector. 

For i = 1, 2 we define two operators G?(z) : L 2 (IR) — > L 2 whose integral kernels are given by 

Gt(z){s,u;s') :=R £ (z)(s,u-s',d) 

G £ 2 {z)(s,u;s') := R £ (z)(s,u; s' , -d) 
and Q £ {z) : L 2 (dVL') -> L 2 given by 

(4.12) Q £ (z)q : 

The operators G\{z + /i n /<5 £2 ), and therefore also Q £ {z + fi n /5 e2 ), are uniformly bounded in e 
for all n = 0, 1, 2, . . . 

(4.13) 





Q £ [z + 



F 2 



Q 



L 2 



the proof of this statement is in appendix |A| 
We introduce also the operators G £ j(z) : L 2 ( 
are given by 



^ c lkl|L2(cK7') , 



L 2 (IR) for i,j = 1,2 whose integral kernels 



Gl :1 (z)(s;s') :=R £ (z)(s,d;s',d) 
Gl, 2 (z){s;s') :=R £ (z)(s,d- s',-d) 
G £ 2>1 (z){s;s') :=R £ (z)(s,-d;s',d) 
G £ (z)(s;s') :=R^z)(s,-d; s',-d) 
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and the operator r e (2;) : L 2 (dQ') — > L 2 {dVl') given by 

'GUz) GUz) 



(4.14) 



The operators G\Az + jj, n /5 e2 ) and, consequently, T £ (z + fi n /6 e2 ), are uniformly bounded in e 



(4.15) 



T £ {z + Vn/S £2 )g\\ L 2 mi) < c||g|| i 2( a n') • 



For the proof of this statement one can refer to the appendix. 

The resolvent R £ (z) can be written in the following form, take ip G L 2 then 

(4.16) R%z)iP = R £ (z)iP + G £ (z)q £ 
where q £ is defined by 

(4.17) q £ = -A £ {z)g £ * (z) A £ {z) := (I + (a- a e )T £ (z))" 1 (a - a £ ) 
where * denotes the adjoint and (a — a £ ) is multiplication operator 

(a — a\) 



[a — a 







(a - a £ 2 ) 



The operator (I+(a— a £ )T £ (z)) 1 in (4.17) is well defined by its Neumann series for e sufficiently 
small, see appendix |A} Now we go back to the proof of (4.9). Using (4.16) we have 



(4.18) 
with 



R £ z + 



5 s2 



F n = Rh[z + 



6 e2 



-A £ z + 



In appendix |A| it is proved that for < e < eq, 



z + 



fie 2 



(4.19) 



\gl\\L 2 (dn') < c~ 



5 e 



We first study the term of second derivative coming from the resolvent Rq{z). The following 
equality holds 



d 2 



■R £ (z; s,u, s',u) 



2 



xRniz: s, u, s', u 



where Rq{z\ s, u, s', u') is the integral kernel of Rq(z). Therefore integrating by parts we have 



du'- 



(s,u) 



ds du 



2 



:Rn(z: S, It, s', u' 



- [R £ (z)\ £ n( p £ n f}(s,u) + J ds'Rl{z-s^s\d){a\{s') - aW n {s\d)f{s') 

+ J ds'R £ (z;s,u,s',-d)(a £ ,(s') - a)(f) £ n (s' , -d)f(s') 

- [Rl(z)X £ n ^f](s,u)+ [Q £ {z){a*-aW{d)f]\ 
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where we used the notation 4> £ (s, d) = ((j) n (s, d), 0^(s, — d)). Then we have 

1 
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d 2 



ds- 



i d 2 



f 



5 s2 dv? e 2 
1 1 2 (-M 



4 



4 



z)R £ (z + l x n /S £2 )cj> £ n f 



+ Rl(z + fi n /6 £2 ) 



: 2^ (K, ~ ^n) 
fie 2 



J-G £ (z + iin/5 £2 )(a £ -aW(d)f 



Since H'yllz, 00 (k) < c an d a > 1, then there exists £o > such that, for all < e < Eq, 
\\v\\l°°{m) ^ cE a ~ l and therefore the perturbative expansion (3.8) can be applied. This implies 
that 



5 s2 



^ c 



.2 ' 



Moreover from ||(a e — a)||i<x>(sn') ^ c5 e /e and estimate (4.13) we have that there exists e such 
that for all < e < e 



(4.20) 



9 s 



5^2 



^ c 1 



L 2 



1 5 £ 

"2 + _ 



^ c- 



L 2 ( 



Now we estimate the second term coming from the r.h.s. of (4.18). We can write it in the 
following equivalent way 

Hn 

§e2j 2* 

We introduce the notation 

r £ d 2 ^s/e) 2 
ds 2 4e 2 ' 

Let us notice that h £ commutes with Rq(z) and T £ (z). Using the commutation property of h e 
we have 



(4.21) 



fie 2 



If 



By estimate (4.13) and (4.19) one can see that the first term on the r.h.s. of (4.21) is bounded 
by c5 e /£ 3 ||/||i2( K ). The second term of (4.21) can be written in the following way 



Q £ [z + 



Q £ lz + 



5 e2 

fin 
fie 2 



h £ - zj A £ \Z + 
h £ - z) A £ (z + 



F 2 

fl n 
fie 2 



Z + 



h £ -z 



fin 

S e2 



If 



z)g e * z + 



fl n 
fie 2 



By an argument similar to the one used in the proof of estimate (4.20) one can prove that 

1 



(4.22) 



h £ -z\g £ * z + 



fin 
6 e2 



L 2 (dn') 



< c- 



Moreover the following identity holds 

h £ -z) A £ (z+j£\ (h £ -z~ X 
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Using the Leibniz rule, we see that 



h £ - z fa - a £ ) Iff - z 



-i 



^ c 



Then there exists Eq such that for < s < eq 



5 £ 5 £ 5 £ 
e e z e 4 



^ c- 



(4.23) 



h £ -z)A £ (z + 



6 e2 



If 



-1 



^ c- 



Therefore it follows from estimates (4.13) and (4.22) and (4.23) that 

d 2 



(4.24) 



Q £ [z + 



F 2 



Q 



^ c- 



%(L 2 (dn'),L 2 (dn')) 
d 
5 £ 



3 ' 



L 2 



Then for a > 3, (4.9) follows from (4.24) and (4.20). Interpolating (4.9) and using the L z 
boundedness of (H £ — Jfj — z)~ l f<fi £ n we immediately obtain 



(4.25) 



d 



ds\ Re fc 2 z 



^ c- 



L 2 ( 



L 2 



The proof of flO| ) comes from p~7j ), ( [48] ), ( |4~9~| and ( |4.25| ). 



□ 



,71 = 0,1,2,... 



Let us consider the family of self-adjoint operators /i^ : if 2 (IR) C L 2 (IR) — > L 2 
defined in (4.1 ). The following lemma concludes the first step in the proof of theorem [TJ it shows 
that in the limit only the diagonal elements of the reduced resolvent survive and that the leading 
term in the reduced Hamiltonian is hi. 



Lemma 2. Let 7 G C, 



00 1 
\ 



) and a > 3 then for all m, n — 0, 1, 2, . . . and for all z G C\M ; 
u - lim (f £ {z) - 8 m ,n{K - z)' 1 ) = . 



Proof. It is sufficient to prove that for all m, n = 0, 1, 2, . . . , /, # G C^°(R) and z G C\M, there 
exists £0 > such that, for all < e < Eq, 



(9, (r E m , n {z) ~ 5 m ,n( h n ~ Z ) *)/) 



L 2 ( 



^- „ _a— 3 1 1 _ 1 1 

^ \\g\\L 2 ( 



L 2 ( 



It is convenient to introduce an intermediate Hamiltonian h n which is the compression of H £ 
to the subspace generated by (p £ n . We define 



hi = -- 



1 ^{s/e) , X £ n {s)-^ n 



+ 



5 



e2 



Since 
(4.26) 



I (#5 (f'm,n( Z ) $m,n(h n Z) )f) L 2m\\ 
<\(9, (r £ m ,n( Z ) ~ $m,n(K ~ ^^f) L 2 (R) | + <W| (g, (K ~ Z )^ ~ (K ~ z )~ X )f) L 2, 



it is sufficient to estimate separately the two terms on the right hand side of (4.26). 
We notice that 

0n(a £ (s)) =£ n + <p £ n ( s ) 

where £ n are the eigenfunctions of the symmetric Robin Laplacian in (—d,d), see (3.6) - (3.7) 
and tp £ n {s) by elementary calculus is given by the following line integral in IR 2 



(4.27) 



V^nfV) • dy_ 



a,a £ (s)] 
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where with a small abuse of notation we have denoted the segment in M 2 between a and a £ (s) 
by [a, a £ (s)}. Notice that 



(4.28) 



\o6 - a || L» (k) < c5 £ je 



which implies Sa = sup sgM \a — a e (s)| ^ c5 £ /e. The eigenstates (j) n {v) are L 2 ((— d, d))- valued 
smooth functions of v G M 2 , in particular \\4> n (v)\\L 2 ((-d,d)) = 1 an d for all «6l 2 and < e < Sq 
there exists an open ball B(a, R £ ) cl 2 with center in a and radius R £ < c5 £ je such that 

(4.29) 



sup \\d Vi (j) n (v)\\ L 2 {{ _ d4)) <^ c , sup \\d Vi d V] <p n (v)\\ L 2 {{ _ d4)) ^ c 



for i,j = 1,2. 

Let T be a given bounded operator in L 2 and let T(s,u; s',u') be its integral kernel. For fixed 
/ G L 2 (IR) we introduce 



C 2 / 0> M ) 



ds'dv! T(s, u\ s', u')cf) n (a £ (s'), u')f(s') 
ds'du' T(s, u; s' , u')(f> n (a £ (s) , v!) f (s') 



and we want to prove 



IICf-CllU^^Vell/IlL^ 



It is convenient to introduce also 



ds'du T(s, u; s', u')£ n (u')f(s') 



and separately estimate ||C/ — Co IU 2 an d IIC2 — Co IU 2 - The estimate of ||Cf — Co Hl 2 trivially 
follows from the boundedness of T and (4.29): 

HCf-Qflb < \\T\\aW)\K<P n -Sn)f\\v <c6°/e\\f\\v m . 
The estimate of ||C/ — Co IU 2 requires a more careful analysis: interchanging integrals, we have 

(4.30) ^C 2 ~~ Co j ( s ' u ) = I dv \ ds'du T(s,u; s' \u')d v <j) n (v, a,u')f(s') + 



+ 



+ 



dv I ds'du T(s,u; s' ,u')d v (j) n {oc, v, u')f(s') + 



a£(s) 



dv\ I dv2 \ ds'du T(s,u; s' ,u')d Ul d l/2 (j) n (i'i,i'2,u')f(s') 



Using (4.28), (4.29) and (4.30) it is straightforward to prove 



IICI-Co^U^c^IItH^l^II/IIl^ 



Let us come back to the proof of (4.26). We notice the following useful identity 



fim,n{9,(h n z) /) L 2 (R) 

with T = ( H £ — fi n /5 e2 — z^j . Therefore we have 
(9, {r £ min (z) - S mjn (h £ n - z)' l )f) L2 



t 1 



■'mi S2 



L 2 



9<Pm, (([ - C 2 ') 



L 2 
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and from the Cauchy-Schwarz inequality we get 



(4.31) 



(9, (r £ m ,n( Z ) ~ 5 m,n(K ~ *) *)/) 



L 2 



^lblU 2 (R)IICi - C 2 / ll < c6 £ /e\ Im A 1 \\g\ 



and this concludes the estimate of the first term at the right hand side of (4.26). 

Let us consider now the second term at the right hand side of equation (4.26). From the first 

resolvent identity we get 



g^ik-zr'-iK-z)- 1 )!) 



LH 



( 2 )^ ^ 



& 



c 2 



By the Cauchy-Schwarz inequality we have 
(4.32) 



IMCfc-*)- 1 -^-*)- 1 )/) 



C I Im z\ 



A £ n (-)-/^n-AlV(-) 2 



C 2 



Since ||7||l°°(r) < c and a > 1, then there exists e$ > such that, for all < e < Eq, 
\\v\\l°°(r) ^ cs a ~ l and therefore the perturbative expansion (3.8) can be applied. This implies 
that 

^(_\2\\ _ ._. ^ ^3(a-l) 



Then we have 



iiA^(-)-^-Ai 2 V(-riu^.^^ 

A £ „(-)-^n-A( 2 V(-) 2 



S e2 



^ CE 



a— 3 



From the last estimate and from equation ( ]4.32 ) we get 

\(g,((K- 



^£ a - 3 c|Im(^)|- 2 ||/|| L2(M) ||^|| L2( 



that together with equation (4.31) concludes the proof of the lemma. 
Proof of theorem [T] 

The proof of theorem [I] follows directly from lemma [TJ lemma [2] and proposition ([2]). 



□ 
□ 



5. Small deformations of the curvature 

In this section we give a generalization of the previous results: we show that deformations of 
the order of e of the angle 9 can lead to a more general coupling in the vertex. 
Let us consider the following scaling for the curvature 

(5.1) 7« - ^ (J) = ^f^7® £ >0, 

where b is a real constant. With this scaling the angle 9 between the straight parts of the curve 
h and I2 is 

-2\ 



9 e = f T{s)ds = Vl + 2Eb9 = (1 + Eb)9 + 0(e 2 



then the scaling (5.1 ) can be interpreted as a deformation, of order e, of the geometric parameter 
9. 

Consider the family of one dimensional Hamiltonians with scaled potential of the form 
(5-2) £^_|l + I + ^ (s/e) . 



GRAPH-LIKE MODELS FOR THIN WAVEGUIDES WITH ROBIN BOUNDARY CONDITIONS 21 

If exists a zero energy resonance for the Hamiltonian h = + v one can define two constants 
c_ and C_|_ cLS it was done in (13 . 11 h . The family of Hamiltonians 

®(hr) ■= {/ G H 2 (R\ 0) s.t. c„/(0 + ) = c+/(0-) , 

c + f (0+) - c_/'(0-) = S(c_/(0") + c + /(0 + ))} 

where 

b:=b f v(s)(f r (s)) 2 ds 

and 

(5.4) ~ hr f:=-^ s^O, 

individuates a family of self-adjoint perturbations of the Laplacian in dimension one (see, e.g., 
0)- 

The following proposition generalizes the result stated in proposition [2j the proof can be read 
in [7], theorem 3.1. (see also [3]). 



Proposition 3. Take h £ and h defined as above and assume that v satisfies conditions (3.10). 
Then two cases can occur: 

(1) There does not exist a zero energy resonance for the Hamiltonian h, then 

u-\im(h £ - z)- 1 = (h - z)- 1 zeC\R. 

(2) There exists a zero energy resonance for the Hamiltonian h, then 

u-\im(h £ - z)- 1 = (h r - z)- 1 zEC\R. 



Let us denote by Q £ the family of domains obtained from Q by scaling 7 and d as stated in (5.1 ) 



and (2.8). Following what was done in section |2| one can define the Robin Laplacian on the 
family of domains Q £ , the operator — A~ e is unitarily equivalent to the operator H £ obtained 
via the substitution 7 — > j £ = ^/l + 2eb^ in (2.9). 

d 1 d 1 d 2 1 ~ 

TTE ._ |_ _Y £ (v 1l) 

ds{l + u^ £ {s)) 2 ds 5 £2 du 2 e 2 K ' h 



with 



~ _ ^ £ (s/e) 2 6 £ /eur"(s/e) 5 (5 £ /e) 2 u 2 r' (s/ef 

V \S,U) .— . ~ r f„\\9 



4(1 + U7] £ (s)) 2 2(1 +U7] £ (s)) 3 4 (l + urj £ (s)) 

and with domain 



T 2fcM\ - , Q ^ < - -A , ^El.\.,.l. -A « 



@(H £ ) := G H\Q!) s.t. -Ms, d) + Si(s)V>(s, d)=0,-^-(s, -d) + a e 2 (sU(s, -d) = 
[ ou ou 

here c?f(s) and 5|(s) are given by 

(5.5) al(s) := a - —^jl—r , a £ 2 (s) := a + ^ 



2(1 +drj E (s)) ' zy J ' 2(1 -dtj £ (s)) 

and we have introduced r/ e (s) := 5 £ / e~f £ (s / e) . 

It is easy to understand that a slightly different version of lemma [T] and lemma [2] holds for the 
Hamiltonian H £ . 
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With the scaling (5.1) the main contribution to the longitudinal part of the Hamiltonian, once 



the dynamics has been reduced to the n-th transverse mode, is given by 
that via proposition [3] leads to a slightly different version of theorem [I] 

For all z G C\M and for all m, n = 0, 1, 2, . . . , let us denote by r e m n (z) the reduced resolvent of 
H e , i.e., the operator in ^(L 2 (IR), L 2 (K)) defined via its integral kernel by 

Ki^ 3 , 8 ') J J 0mO, M ) (h e ~ ~ z ) (s,u;s',u')(f) £ n (s',u)dudu . 

Theorem 2. Assume that T has no self-intersections and that 7 G C^°(R) ; moreover take a > 3 
then for all n, m = 0, 1, 2, . . . two cases can occur: 

(1) For all n such that there does not exist a zero energy resonance for h n we have 

u- limr^ n (z) = S m Jh - z)' 1 z G C\R, m, n = 0, 1, 2, . . . . 

(2) For all n such that there exists a zero energy resonance, f r ^ n , for h n we have 

u-limr^ n O) = 5 min (h rtn - z)' 1 z G C\R, m, n = 0, 1, 2, . . . . 

e— >0 



where h rjU is defined according to (5.3) and (5.4). 



6. Conclusions and remarks 

We have studied the convergence of the Robin Laplacian on a waveguide to an operator on a 
graph made up of two edges and one vertex. In this setting we were able to give detailed results 
on the convergence in norm resolvent sense. Our analysis takes into account the projection on 
each transverse mode and distinguishes generic and non-generic cases leading respectively to 
decoupling and non-decoupling gluing conditions in the vertex. Non-generic cases are related 
to the existence of zero energy resonances for the leading term of the effective Hamiltonian. 
The relation between gluing conditions in the vertex and the resonance is rigorously stated in 



equations (3.11) - (3.13) 



The existence of a zero energy resonance is an exceptional event and in general it is destroyed 
by slight deformations of the potential. As it is clearly shown in figure [2j apart for few special 
cases, the coefficients (3 n change as n changes. For these reasons even in the simpler example 
of a waveguide in most of the cases the operator on the corresponding graph is defined by 
conditions in the vertex of decoupling type and in general, for fixed a, one can expect coupling 
at most in one transverse mode. This is in agreement with previous results derived in [3], [7], 
[30] and [15]. 

It is interesting to note that in this simple model all the interplay between geometry and 
boundary conditions on the initial domain, is reduced to the value of (3 n . In particular a positive 
sign of p n , may rule out the possibility of having a zero energy resonance giving decoupling 
conditions in the limit. We conjecture that much more complicate geometries than the vertex 
region of a strip with constant width may give different threshold singularities opening up for 
the possibility of much more general non decoupling conditions. 
We discuss with more detail our result for some suitable choices of the parameter a. 
Formally the limit a ±00 gives Dirichlet conditions on the boundary of the waveguide. 
This case was discussed in two former works [3 J and [7], the results stated there are included 



in our paper. This can be seen by studying the asymptotic behavior of equation (3.15) and 



in particular by noticing that (3 n (±oo) = —1/4, see also figure^! The anomalous behavior of 
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Po(a) and (5\{a) for a — > — oo is related to the existence of two negative eigenvalues for h a when 
a < —1, see the plot of /U (a) and Hi(a) in figure [I] It is important to notice that Dirichlet 
boundary conditions on dfl are preserved by the unitary map (2.6). As a consequence the 
analysis of the transverse part of the Hamiltonian is particularly simple in this case: transverse 
modes and eigenvalues do not depend on s and the Hamiltonian on the graph does not depend 
on n. 

The case a = reproduces Neumann boundary conditions. As expected even in our model 
the Neumann waveguide, on the ground state, is approximated by free conditions in the vertex 
and this behavior is summarized by /3q(0) = 0. The reduced Hamiltonian relative to the lowest 
transverse mode in the Neumann case is the free Laplacian on the line and the zero energy 
resonance is the constant function. In this case the constants c_ and c + coincide and the limit 
Hamiltonian is defined by free conditions in the vertex. 

Since in the Neumann case /3 n (0) = 3/4 for n — 1,2, ... , the Hamiltonian on the graph obtained 
by projecting onto the excited transverse modes has decoupling conditions in the vertex. 
We want to stress an important difference between our approach and the one generically used 
in other works. In our setting the domain Q e is defined via two characteristic lengths: the 
width of the waveguide 5 e and the range e on which the curvature varies. Since we assume that 
for small e we have 5 £ <C e, we can make use of the adiabatic separation of the dynamics that 
leads in a natural way toward an analysis of the problem in two steps. The existence of two 
different scales of length makes it difficult to compare our model with the one generically used 
in the works cited in the introduction in which the scaling of the network in a neighborhood of 
the vertex is assumed to be isotropic. 

The analysis of the case 5 s = e would be of great interest because it reproduces the scaling of 
[To] and [SO]. The failing of the adiabatic approach makes this case much more complicated 
and we guess that in this setting the uniform resolvent convergence could be too demanding. 
Even in more complicated settings, such as graphs with three or more edges, it seems reasonable 
to expect that deformations of order e of networks associated to non decoupling conditions in 
the vertex could lead to more general gluing conditions. This idea is suggested from the result 
stated in section [5] and was already envisaged in [7j. 
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Appendix A. Proof of formula (4.16) 



We mark that formula (4.16 ) can be proved by making use of a very general technique developed 
by A. Posilicano in [31 J and [32] • For convenience of the reader we give a more direct proof 
of the formula in this appendix together with the proof of several technical estimates used in 
lemma [TJ 

Let us start with the proof of formula (4.16). We denote by Q £ the quadratic form associated 
to the operator Hq defined in (4.10) - (4.11), it is given by 



Q £ o[<p,i>} 



17 2 _, 
72 — ^ 



ds du 



+ a 



dip dtp 1 dcp dip 
ds ds 5 s2 du du e 2 4 

(jp(s, d)ip(s, d) + Tp(s, —d)ip(s, —d))ds . 
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We recall also that the quadratic form Q e associated to H £ defined in (4.2) - (4.3) is 



(A.l) Q £ [y^] = Q E [<p,i(>] + / [(al(s)-a)ip(s,d)^(s,d) + (al(s)-a)(p(s,-d)ij(s,-d)}ds^ 

Jm. 

Let us prove that for all i/j G L 2 

R £ (z)i; = R £ (z)^P + Q £ {z)q £ 



where Q £ {z) was defined in (4.12) and q £ is a solution of equation 
(A.2) q £ + (or - of) {T £ {z)q £ + G £ *{z)^) = . 

For ip G L 2 , we define 

T £ (z)^ = R £ (z)^ + g £ (z)q £ 
and we look for conditions on q e such that 

(A.3) Q*[<p, T £ (z)^} - z(ip, T £ (z)i;) = {cp, i>) 



for any <p G C£°(IR 2 ). If ( fO) ) holds then T e (z) = R £ (z). If we start from ( |Al| , integrate by 
parts and notice that Ro(z)if) G @(Hq), we have 

Q £ [if,T £ (z)^-z(if,T £ (z)^) = + [f + («-« £ ) (rW + S'W)])^*,,) 

with <£>(s) = (<£>(s, d), cp(s, —d)), which implies (A.2). Now we have to prove that equation ( A.2[ ) 
has a unique solution for all ip G L 2 . To this aim we need to prove that T e (z) and G £ {z) are 
bounded. Let us prove estimates (4.13) and (4.15) that are required in the proof of lemma [Tj 
then the boundedness of T £ (z) and Q £ {z) will be an obvious consequence. 
To prove ( |4.13 ) we need to prove that that G\ (z + /i n /5 e2 ) are operators uniformly bounded in 
e for all n = 0, 1, 2, ... ; we shall prove the boundedness only of G\{z + fi n /5 e2 ), the proof for 
G% { z + Jfj) is similar and will be omitted. 

d 2 1 7 2 (s/e) 

We denote by V £ (d\) the projector valued measure associated to the operator — -— — 

CIS S Q 

and by a £ its spectrum. Given a function q G L 2 (M), by the spectral theorem we have 



L2 J^j^ \ X + (Vj-Vn)/5 e2 -z\ 

It is easy to see that for Im z ^ 



:(q,V £ (d\)q) LH 



this is proved by decomposing the series into two parts 

oo c2(J\ c2i 



« c, 



3=0 



e 3 (d) 



A + (fij -fi n )/5 £2 



- Z\ 2 ^ ^ E Im _ 



| Mi - Mn + ^ 2 (A-z)p 



^| Mi - Mn + ^( A -^)|2 



The first term in the right hand side can be easily estimated by the imaginary part of z 



* e4 £ 



\fij — jji n + 5 e2 (X — z)\ 2 I Im z\ 



To estimate the last term we notice that, since A G a £ C [— c/e 2 , oo) there exists £o such that 
for all < e < Eo and for all j > n 

\fij ~ Mn + 5 e2 (A - z)\ ^ - //„ - c(<f /e) 2 | 
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from which it follows that for all < e < Sq and uniformly in A G o £ 

8(d) 



E 

j>n 



\fij ~ H n + 6 e2 (\- z)\' 



^ c. 



Therefore from the definition of spectral projection we get 



Q 



L 2 



< c \\q\\lH 



and then G\ (z + is bounded uniformly in e. The boundedness of Q £ {z) comes from a 
similar argument and the proof will be omitted. The same estimates obviously hold for Q £ *{z). 
By noticing that for all q G L 2 (M) 

2 



G 1A Z + 



5 £2 , 





2 












L 2 m 





E 

j=0 



X + (flj - H n )/& 



e2 



L 2 ( 



dX. 



and by an argument similar to the one used before one can show that for all Im z ^ there 
exists an e such that for all < e < e 

2 

< c. 



oo 



z 



From which it follows that G\ x (z + fi n /5 e2 ) is uniformly bounded in e. The same holds for all 
Gfj(z + fi n /8 £2 ) and, consequently, for T £ (z + fi n /5 e2 ). The boundedness of T £ (z) easily comes. 
Now we can go back to formula (A. 2 ) and prove that there exists a unique solution q £ G L 2 (dQ'). 
We have just proved that there exists e such that, for < e < e , \\(a — a e )r e (z)|| L 2( 9n /) < 1, 
then (J + (a — a £ )T £ (z))" 1 is well defined by its Neumann series and 

q e = -(I + (a - a £ )T £ (z)y\a - a £ )g £ *(z)^j 



is the solution of (A.2). 

This concludes the proof of formula 4.16 From the inequality \\(a — cr^lli^dn') ^ 3 e /z we also 
get lk £ ||L 2 (9n) ^ c^/^II^IIl 2 and we proved all the estimates that we used in lemma [TJ 

References 

[1] R. R. Adams, Sobolev spaces, Academic Press, New York, 1975. 

[2] S. Albeverio, Z. Brzezniak, and L. Dabrowski, Fundamental solutions of the heat and Schrddinger equations 

with point interactions, J. Funct. Anal. 130 (1995), 220-254. 
[3] S. Albeverio, C. Cacciapuoti, and D. Finco, Coupling in the singular limit of thin quantum waveguides, J. 

Math. Phys. 48 (2007), 032103. 
[4] G. Berkolaiko, R. Carlson, S. Fulling, and P. Kuchment, Quantum graphs and their applications, Contem- 
porary Math., vol. 415, American Math. Society, Providence, R.I., 2006. 
[5] D. Bolle, F. Gesztesy, and S. F. J. Wilk, A complete treatment of low- energy scattering in one dimension, 

J. Operator Theory 13 (1985), no. 1, 3-32. 
[6] G. Bouchittc, M. L. Mascarenhas, and L. Trabucho, On the curvature and torsion effects in one dimensional 

waveguides, ESAIM Control Optim. Calc. Var. 13 (2007), no. 4, 793-808. 
[7] C. Cacciapuoti and P. Exner, Nontrivial edge coupling from a Dirichlet network squeezing: the case of a 

bent waveguide, J. Phys. A: Math. Theor. 40 (2007), no. 26, F511-F523. 
[8] G. Dell' Antonio and L. Tcnuta, Quantum graphs as holonomic constraints, J. Math. Phys. 47 (2006), 

072102. 

[9] J. Dittrich and J. Kffz, Curved planar quantum wires with Dirichlet and Neumann boundary conditions, J. 
Phys. A 35 (2002), no. 20, L269-L275. 
[10] P. Duclos and P. Exner, Curvature- induced bound states in quantum waveguides in two and three dimen- 
sions, Rev. Math. Phys. 7 (1995), no. 1, 73-102. 



26 CLAUDIO CACCIAPUOTI AND DOMENICO FINCO 

[11] P. Exner and O. Post, Convergence of spectra of graph-like thin manifolds, J. Geom. Phys. 54 (2005), 
77-115. 

[12] P. Exner and P. Seba, Free quantum motion on a branching graphs, Rep. Math. Phys. 28 (1989), no. 1, 
7-26. 

[13] M. I. Freidlin and A. D. Wentzel, Diffusion processes on graphs and averaging principle, Ann. Probab. 21 
(1993), no. 4, 2215-2245. 

[14] P. Freitas and D. Krejcink, Location of the nodal set for thin curved tubes, Indiana Univ. Math. J. 57 
(2008), no. 1, 343-376. 

[15] D. Grieser, Spectra of graph neighborhoods and scattering, arXiv:0710.3405v4 [math.SP] (2007). 

[16] M. Harmer, Hermitian symplectic geometry and the factorization of the scattering matrix on graphs, J. 

Phys. A 33 (2000), no. 49, 9015-9032. 
[17] P. Hejci'k and T. Chcon, Irregular dynamics in a solvable one- dimensional quantum graph, Phys. Lett. A 

356 (2006), no. 4-5, 290-293. 
[18] N. E. Hurt, Mathematical physics of quantum wires and devices., Kluwer Academic Publishers, Dordrecht, 

2000. 

[19] A. Jensen and G. Nenciu, A unified approach to resolvent expansions at thresholds, Rev. Math. Phys. 13 
(2001), no. 6, 717-754. 

[20] T. Kato, Perturbation theory for linear operators, Springer- Ver lag, Berlin Heidelberg New York, 1980. 
[21] V. Kostrykin and R. Schrader, Kirchhoff's rule for quantum wires, J. Phys. A: Math. Gen. 32 (1999), no. 4, 
595-630. 

[22] D. Krejcifik, Spectrum of the laplacian in a narrow curved strip with combined dirichlet and neumann 
boundary conditions, ESAIM: Control, Optimisation and Calculus of Variations (2008), to appear. Preprint 
available on |arXiv:0803.0937| [math.SP]. 

[23] D. Krejcink and J. Knz, On the spectrum of curved quantum waveguides, RIMS, Kyoto University 41 
(2005), no. 3, 757-791. 

[24] P. Kuchment, Graph models for waves in thin structures, Waves Random Media 12 (2002), no. 4, R1-R24. 

[25] , Quantum graphs. I. Some basic structures, Waves Random Media 14 (2004), no. 1, S107-S128. 

[26] , Quantum graphs. II. Some spectral properties of quantum and combinatorial graphs, J. Phys. A: 

Math. Gen. 38 (2005), no. 22, 4887-4900. 

[27] P. Kuchment and H. Zeng, Convergence of spectra of mesoscopic systems collapsing onto a graph, J. Math. 
Anal. Appl. 258 (2001), no. 2, 671-700. 

[28] K. Kuwae and T. Shioya, Convergence of spectral structures: a functional analytic theory and its applica- 
tions to spectral geometry, Comm. Anal. Geom. 11 (2003), no. 4, 599-673. 

[29] J. T. Londergan, J. P. Carini, and D. P. Murdock, Binding and scattering in two-dimensional systems, 
vol. 60, Springer- Ver lag, Berlin, 1999, Lecture Notes in Physics Monographs. 

[30] S. Molchanov and B. Vainberg, Laplace operator in networks of thin fibers: spectrum near the threshold, 
|arXiv:0704.2"795y i [math-ph] (2007). 

[31] A. Posilicano, A Krein-like formula for singular perturbations of self-adjoint operators and applications, J. 
Funct. Anal. 183 (2001), 109-147. 

[32] , Self-adjoint extensions of restrictions, arXiv:math-ph/0703078v2 (2007). 

[33] O. Post, Branched quantum wave guides with Dirichlet boundary conditions: the decoupling case, J. Phys. 

A: Math. Gen. 38 (2005), no. 22, 4917-4931. 

[34] , Spectral convergence of quasi-one-dimensional spaces, Ann. Henri Poincare 7 (2006), 933-973. 

[35] M. Reed and B. Simon, Methods of modern mathematical physics, vol. II: Fourier Analysis, Self-Adjointness, 

Academic Press, New York-San Francisco-London, 1975. 
[36] J. Rubinstein and M. Schatzman, Variational problems on multiply connected thin strips. I. Basic estimates 

and convergence of the Laplacian spectrum, Arch. Ration. Mech. Anal. 160 (2001), no. 4, 271-308. 
[37] K. Ruedenberg and C. W. Scherr, Free-electron network model for conjugated systems. I. Theory, J. Chem. 

Phys. 21 (1953), no. 9, 1565-1581. 
[38] Y. Saito, Convergence of the Neumann Laplacian on shrinking domains, Analysis (Munich) 21 (2001), 

no. 2, 171-204. 

[39] S. Sclbcrhcrr, Analysis and simulation of semiconductor devices, Springer- Ver lag, Wien New York, 1984. 



GRAPH-LIKE MODELS FOR THIN WAVEGUIDES WITH ROBIN BOUNDARY CONDITIONS 

Cacciapuoti: Czech Technical University, Doppler Institute 
Current address: Bfchova 7, 11000 Prague, Czech Republic 
E-mail address: cacciapuoti@ujf.cas.cz 

Finco: Department of Mathematics "G. Castelnuovo" , University of Rome "La Sapienza" 
Current address: P.le A.Moro 2, 00185 Rome, Italy 
E-mail address: finco@mat.uniromal.it 



